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DETAILED ACTION 

1 . The response filed 6/1 0/1 0 has been received and considered. Claims 1-17, 1 9- 
29, 31-33, 37 and 38 have been presented for examination. Claims 37 and 38 are new. 

Claim Rejections - 35 USC § 102 

2. The following is a quotation of the appropriate paragraphs of 35 U.S.C. 1 02 that 
form the basis for the rejections under this section made in this Office action: 

A person shall be entitled to a patent unless - 

(e) the invention was described in (1 ) an application for patent, published under section 1 22(b), by 
another filed in the United States before the invention by the applicant for patent or (2) a patent 
granted on an application for patent by another filed in the United States before the invention by the 
applicant for patent, except that an international application filed under the treaty defined in section 
351(a) shall have the effects for purposes of this subsection of an application filed in the United States 
only if the international application designated the United States and was published under Article 21(2) 
of such treaty in the English language. 

3. Claims 1-7, 12-15, 17, 19-22, 25-29, 31, 32, 37 and 38 are rejected under 35 
U.S.C. 102(e) as being anticipated by Shah et al (US Patent Application Publication 
2004/0084180). 

4. As to Claim 1 , Shah et al teaches: a method of determining production rates in a 
well (paragraph 0005, "...methods and systems for estimating multi-phase fluid flow 
rates in a subterranean well..."; paragraph 0009, "...estimates multi-fluid flow rates are 
provided for a plurality of selected well locations"), comprising: determining a model of 
temperature as a function of zonal flow rates in the well (paragraph 0025, lines 4-6, "It is 
well known that each of the entries 24 has its own fluid phase (whether it be oil, water or 
gas), flow rate, temperature and hydrocarbon mixture composition"; paragraph 0027, 
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lines 1-14, "...the model 30 takes into account the conservation of energy and mass, 
and consequently simulates the evolution of the temperature of the flowing fluid... It is 
known that pressure and temperature of the flowing fluids change as they travel up or 
down a flow path... Therefore, model 30 should take such factors into consideration"; 
paragraph 0029, lines 8-11, "The user of the invention may specify... the temperature at 
each entry point 24 within the well 10"; paragraph 0036, lines 9-12, "In step 106, the 
mathematical model for the wellbore 12 is run to calculate the expected pressure and 
temperature values at the downhole sensor locations..."); measuring temperatures at a 
plurality of locations in the well (paragraph 0015; paragraph 0025, lines 16-24, "...a 
plurality of sensors 27 are deployed to take measurements at the various production 
zones 26 or other production points of interest inside the wellbore 12..."); and inverting, 
by a computer, the measured temperatures by applying the model to determine an 
allocation of production rates from different producing zones in the well (Figure 2 and 
description, paragraph 0030; paragraph 0039, lines 4-5; paragraph 0031, lines 1-4, "The 
data path 31 supplies transient data to the model, 30, such as... temperature data 43 
measured at multiple downhole locations"; paragraph 0036, "Transient well data is 
measured. ..including pressure and temperature data in the wellbore 12 above each flow 
entry being produced... measurement above each flow entry is not required for the 
solution of the inverse problem... In step 106, the mathematical model of the wellbore 12 
is run to calculate the expected pressure and temperature values at the downhole 
sensor locations, the expected volumetric phase flow rates at the wellhead 22, and 
sensitivity coefficients of the model response to each phase flow rate at each fluid entry 
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location"; paragraph 0037, lines 4-10, "...the model-calculated expected pressure and 
temperature values for various well locations of step 106 are preferably compared with 
the measured temperature and pressure values obtained in step 104 with respect to 
those same well locations..."; paragraph 0038, "...the deviation between the calculated 
and measured quantities (of step 108) may be used with the sensitivity coefficients of 
the model (from step 106) to determine changes necessary... In this way, the modeling 
comparisons may be reiterated... until an approximate match (within acceptable 
tolerances) is obtained between the calculated well properties and related flow rates 
and the measured well properties and flow rates.. .if the measured volumetric phase 
rates and pressure and temperature readings are in tolerable agreement with the 
expected values predicted by the model 30, the final estimates of the multiphase flow 
rates are provided as shown in step 118"; claim 27, "...the multiphase fluid flow rate 
program further comprises a model inversion algorithm"; paragraph 0009, "...estimated 
multi-fluid flow rates are provided for a plurality of selected well locations"); wherein the 
inverting comprising using an optimization algorithm that solves an optimization problem 
for calculating the production rates (Figure 3. elements 108, 110; paragraph 0006, lines 
5-8, "The multi-phase fluid flow rates are estimated by iteratively comparing measured 
static and transient well conditions with the model for the well"; paragraphs 0037-0038, 
"...In step 108, the model-calculated expected pressure and temperature values for 
various well locations of step 106 are preferably compared with the measured 
temperature and pressure values obtained in step 104 with respect to those same well 
locations. Thus at step 108, the actual transient data is compared with the calculated 
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expectations of the model...", "...the deviation between the calculated and measured 
quantities (of step 108) may be used with the sensitivity coefficients of the model (from 
step 106) to determine changes necessary. ..In this way, the modeling comparisons may 
be reiterated. ..until an approximate match (within acceptable tolerances) is obtained 
between the calculated well properties and related flow rates and the measured well 
properties and flow rates. ..if the measured volumetric phase rates and pressure and 
temperature readings are in tolerable agreement with the expected values predicted by 
the model 30, the final estimates of the multiphase flow rates are provided as shown in 
step 1 1 8"), where the optimization problem minimizes an error between the measured 
temperatures and corresponding temperatures calculated by the model (Figure 3, 
elements 108, 112; paragraphs 0037 and 0038, "...the deviation between the calculated 
and measured quantities (of step 108) may be used with the sensitivity coefficients of 
the model (from step 106) to determine changes necessary for the estimate of phase 
flow rates at each well entry point. ..the modeling comparisons may be reiterated. ..until 
an approximate match (within acceptable tolerances) is obtained between the 
calculated well properties and related flow rates and the measured well properties and 
flow rates..."). 

5. As to Claim 2. Shah et al teaches: wherein determining the model comprises 
determining the model for a single-phase liquid producing well (paragraph 0026, lines 1- 
7). 
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6. As to Claim 3, Shah et al teaches: wherein determining the model comprises 
determining the model for a multi-layer producing well (Figure 1 ; paragraph 0025, lines 
10-13). 

7. As to Claim 4, Shah et al teaches: wherein determining the model comprises 
determining the model for a multi-layer, single-phase liquid producing well (Figure 1; 
paragraph 0025, lines 10-13; paragraph 0026, lines 1-7). 

8. As to Claim 5, Shah et al teaches: wherein determining the model comprises 
determining the model for a multi-layer, multi-phase liquid producing well (Figure 1; 
paragraph 0025, lines 10-13; paragraph 0026, lines 1-7). 

9. As to Claim 6, Shah et al teaches: wherein measuring the temperatures 
comprises measuring temperature with a distributed temperature sensor (paragraph 
0025, lines 16-24). 

10. As to Claim 7, Shah et al teaches: wherein the inverting comprises determining a 
degree of certainty in the production rates allocated (paragraph 0012; paragraph 0037, 
lines 9-12; paragraph 0038, lines 6-15). 

11. As to Claim 37, Shah et al teaches: measuring a total flow rate of the well at a 
wellhead (paragraph 0011, lines 7-9, "...measuring. ..the wellhead flow rate..."; 
paragraph 0036, lines 8-9, "Volumetric flow rate measurements for each phase at the 
wellhead 22 are also obtained"; Figure 3, step 104; Figure 4, element 70; paragraph 
0041, lines 1-2, "...the wellhead flow rate is also provided to the model..."); and 
allocating, by the model, the total flow rate among the different producing zones based 
on the measured temperatures (paragraph 0031, lines 1-4, "The data path 31 supplies 
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transient data to the model, 30, such as. ..temperature data 43 measured at multiple 
downhole locations"; paragraph 0036, "Transient well data is measured. ..including 
pressure and temperature data in the wellbore 12 above each flow entry being 
produced. ..measurement above each flow entry is not required for the solution of the 
inverse problem... In step 106, the mathematical model of the wellbore 12 is run to 
calculate the expected pressure and temperature values at the downhole sensor 
locations, the expected volumetric phase flow rates at the wellhead 22, and sensitivity 
coefficients of the model response to each phase flow rate at each fluid entry location"; 
paragraph 0037, lines 4-10, "...the model-calculated expected pressure and 
temperature values for various well locations of step 106 are preferably compared with 
the measured temperature and pressure values obtained in step 104 with respect to 
those same well locations..."; paragraph 0038, "...the deviation between the calculated 
and measured quantities (of step 108) may be used with the sensitivity coefficients of 
the model (from step 106) to determine changes necessary... In this way, the modeling 
comparisons may be reiterated... until an approximate match (within acceptable 
tolerances) is obtained between the calculated well properties and related flow rates 
and the measured well properties and flow rates. ..if the measured volumetric phase 
rates and pressure and temperature readings are in tolerable agreement with the 
expected values predicted by the model 30, the final estimates of the multiphase flow 
rates are provided as shown in step 118"; claim 27, "...the multiphase fluid flow rate 
program further comprises a model inversion algorithm"; paragraph 0009, "...estimated 
multi-fluid flow rates are provided for a plurality of selected well locations"). 
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12. As to Claim 12, Shah et al teaches: a method of determining flow rates in a well 
(paragraph 0005, "...methods and systems for estimating multi-phase fluid flow rates in 
a subterranean well..."; paragraph 0009, "...estimates multi-fluid flow rates are provided 
for a plurality of selected well locations"), comprising: measuring temperatures at a 
plurality of points along the well (paragraph 0015; paragraph 0025, lines 16-24, "...a 
plurality of sensors 27 are deployed to take measurements at the various production 
zones 26 or other production points of interest inside the wellbore 12...") having a 
plurality of well zones and a plurality of liquid phases (Figure 1; paragraph 0025, lines 
10-13; paragraph 0026, lines 1-7); measuring a total flow rate from the well (paragraph 
0011, lines 7-9, "...measuring. ..the wellhead flow rate..."; paragraph 0036, lines 8-9, 
"Volumetric flow rate measurements for each phase at the wellhead 22 are also 
obtained"; Figure 3, step 104; Figure 4, element 70; paragraph 0041, lines 1-2, "...the 
wellhead flow rate is also provided to the model..."); and determining, by a computer, 
flow rates of the plurality of liquid phases through each of the plurality of well zones 
(Abstract, lines 7-10, "Multi-phase fluid flow estimates may be obtained for the various 
liquid and gaseous fluids in the well (10) at multiple well locations (24)"; paragraph 
0009, "...estimated multi-fluid flow rates are provided for a plurality of selected well 
locations"; Figure 2 and description, paragraph 0030; paragraph 0039, lines 4-5; 
paragraph 0036, lines 9-14, "...sensitivity coefficients of the model response to each 
phase flow rate at each fluid entry location..."; paragraph 0037, lines 1-4, "...the 
expected wellhead volumetric flow of each phase calculated in step 106 is compared 
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with the measured volumetric phase flow rate obtained instep 104"; paragraph 0038, 
"...the modeling comparisons may be re-iterated. ..until an approximate match (within 
acceptable tolerances) is obtained between the calculated well properties and related 
flow rates and the measured well properties and flow rates... the final estimates of 
multiphase flow rates are provided...") via the measured temperatures (paragraph 001 1 , 
"...the method of estimating multi-phase fluid flow rates in a subterranean well. ..Steps 
are also provided for measuring the transient temperature, pressure and wellhead flow 
rate. The subterranean well is modeled using these measurements to estimate the 
multi-phase low rates in the subterranean well"; paragraph 0031, lines 1-4, "The data 
path 31 supplies transient data to the model, 30, such as. ..temperature data 43 
measured at multiple downhole locations"; Figure 3, element 104-1 18 and description, 
paragraphs 0036-0038); wherein the determining comprises inverting the measured 
temperatures by applying a model (paragraph 0031, lines 1-4, "The data path 31 
supplies transient data to the model, 30, such as... temperature data 43 measured at 
multiple downhole locations"; paragraph 0036, "Transient well data is 
measured. ..including pressure and temperature data in the wellbore 12 above each flow 
entry being produced... measurement above each flow entry is not required for the 
solution of the inverse problem... In step 106, the mathematical model of the wellbore 12 
is run to calculate the expected pressure and temperature values at the downhole 
sensor locations, the expected volumetric phase flow rates at the wellhead 22, and 
sensitivity coefficients of the model response to each phase flow rate at each fluid entry 
location"; paragraph 0037, lines 4-10, "...the model-calculated expected pressure and 
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temperature values for various well locations of step 1 06 are preferably compared with 
the measured temperature and pressure values obtained in step 104 with respect to 
those same well locations..."; paragraph 0038, "...the deviation between the calculated 
and measured quantities (of step 108) may be used with the sensitivity coefficients of 
the model (from step 106) to determine changes necessary... In this way, the modeling 
comparisons may be reiterated... until an approximate match (within acceptable 
tolerances) is obtained between the calculated well properties and related flow rates 
and the measured well properties and flow rates. ..if the measured volumetric phase 
rates and pressure and temperature readings are in tolerable agreement with the 
expected values predicted by the model 30, the final estimates of the multiphase flow 
rates are provided as shown in step 118"; claim 27, "...the multiphase fluid flow rate 
program further comprises a model inversion algorithm"; paragraph 0009, "...estimated 
multi-fluid flow rates are provided for a plurality of selected well locations"), wherein the 
inverting comprises allocating bv the total flow rate among the plurality of well zones 
(paragraph 0009, "...estimated multi-fluid flow rates are provided for a plurality of 
selected well locations"; Figure 3, element 104, "Volumetric flow rate of each phase at 
the well head"; Figure 4, element 70; paragraph 001 1 , lines 7-1 1 ; paragraph 0036, lines 
8-9, "Volumetric flow rate measurements for each phase at the wellhead 22 are 
obtained"; paragraph 0037, lines 1-4, "...the expected wellhead volumetric flow of each 
phase calculated in step 106 is compared with the measured volumetric phase flow rate 
obtained in step 104"; paragraph 0038, "...the deviation between the calculated and 
measured quantities (of step 108) may be used with the sensitivity coefficients of the 
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model (from step 106) to determine changes necessary for the estimate of phase flow 
rates at each well entry point. ..if the measured volumetric phase rates and pressure and 
temperature readings are in tolerable agreement with the expected values predicted by 
the model 30, the final estimates of the multiphase flow rates are provided as shown at 
step 118"). 

13. As to Claim 13, Shah et al teaches: wherein measuring the temperature at the 
plurality of points comprises utilizing a distributed temperature sensor (paragraph 0025, 
lines 16-24). 

14. As to Claim 14, Shah et al teaches: wherein determining the flow rates 
comprises constructing the model of temperature as a function of zonal flow rates in the 
well (paragraph 0025, lines 4-6, "It is well known that each of the entries 24 has its own 
fluid phase (whether it be oil, water or gas), flow rate, temperature and hydrocarbon 
mixture composition"; paragraph 0027, lines 1-14, "...the model 30 takes into account 
the conservation of energy and mass, and consequently simulates the evolution of the 
temperature of the flowing fluid... It is known that pressure and temperature of the 
flowing fluids change as they travel up or down a flow path... Therefore, model 30 
should take such factors into consideration"; paragraph 0029, lines 8-1 1 , "The user of 
the invention may specify... the temperature at each entry point 24 within the well 10"; 
paragraph 0036, lines 9-12, "In step 106, the mathematical model for the wellbore 12 is 
run to calculate the expected pressure and temperature values at the downhole sensor 
locations..."), and using the model to invert the measured temperatures in allocating the 
flow rates from the plurality of well zones based on the measured total flow rate 
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(paragraph 001 1 , lines 7-9, "wellhead flow rate"; paragraph 0031 , lines 1-4, "The data 
path 31 supplies transient data to the model, 30, such as... temperature data 43 
measured at multiple downhole locations"; paragraph 0036, "Transient well data is 
measured. ..including pressure and temperature data in the wellbore 12 above each flow 
entry being produced... measurement above each flow entry is not required for the 
solution of the inverse problem... Volumetric flow rate measurements for each phase at 
the wellhead 22 are also obtained. In step 106, the mathematical model of the wellbore 
12 is run to calculate the expected pressure and temperature values at the downhole 
sensor locations, the expected volumetric phase flow rates at the wellhead..."; 
paragraph 0037, lines 1-10, "...the expected volumetric flow of each phase calculated in 
step 106 is compared with the measured volumetric phase flow rate obtained in step 
104... the model-calculated expected pressure and temperature values for various well 
locations of step 106 are preferably compared with the measured temperature and 
pressure values obtained in step 104 with respect to those same well locations..."; 
paragraph 0038, "...the deviation between the calculated and measured quantities (of 
step 1 08) may be used with the sensitivity coefficients of the model (from step 1 06) to 
determine changes necessary. ..In this way, the modeling comparisons may be 
reiterated... until an approximate match (within acceptable tolerances) is obtained 
between the calculated well properties and related flow rates and the measured well 
properties and flow rates. ..if the measured volumetric phase rates and pressure and 
temperature readings are in tolerable agreement with the expected values predicted by 
the model 30, the final estimates of the multiphase flow rates are provided as shown in 
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step 118"; claim 27, "...the multiphase fluid flow rate program further comprises a model 
inversion algorithm"; paragraph 0009, "...estimated multi-fluid flow rates are provided for 
a plurality of selected well locations"). 

15. As to Claim 15, Shah et al teaches: wherein determining the flow rates 
comprises determining flow rates of oil and water phases during production (Abstract, 
lines 7-10, "Multi-phase fluid flow estimates may be obtained for the various liquid and 
gaseous fluids in the well (10) at multiple well locations (24)"; paragraph 0009, 
"...estimated multi-fluid flow rates are provided for a plurality of selected well locations"; 
paragraph 0037, lines 1-4, "...the expected wellhead volumetric flow of each phase 
calculated in step 106 is compared with the measured volumetric phase flow rate 
obtained instep 104"; paragraph 0038, lines 11-15, "...if the measured volumetric phase 
rates... are intolerable agreement with the expected values predicted by the model 30, 
the final estimates of the multiphase flow rates are provided..."; paragraph 0035, line 5; 
paragraph 0028, lines 1-4). 

16. As to Claim 17, Shah et al teaches: wherein inverting the measured 
temperatures comprises utilizing an optimization algorithm that solves an optimization 
problem for calculating the flow rates, where the optimization problem minimizes an 
error between the measured temperatures and corresponding temperatures calculated 
by the model (Figure 3 and description, paragraphs 0033-0038). 



1 7. As to Claim 19, Shah et al teaches: a system (Figures 1 , 2), comprising: a 
temperature sensor deplovable with a production completion along a wellbore to sense 
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temperature data at a plurality of wellbore locations during production (paragraph 0025, 
lines 16-24, "...a plurality of sensors 27 are deployed to take measurements at the 
various production zones 26 or other points of interest inside the wellbore 12. The 
sensors 27 are preferably downhole temperature and pressure transducers coupled to 
computer 32. ..may include fiber optic distributed temperature sensing ("DTS") 
systems..."); and a processor system configured to receive the sensed temperature data 
and allocate flow rates from a plurality of wellbore zones based on the sensed 
temperature data (paragraph 0025, lines 16-24, "The sensors 27 are preferably 
downhole temperature and pressure transducers coupled to computer 32..."; paragraph 
0026, lines 5, "...the computer 32 incorporates the functionality of a mathematical model 
30 designed to simulate the physical processes of the flow of multi-phase fluid. ..within 
the wellbore 12"; paragraph 0031, lines 1-4; paragraphs 0036-0038; paragraph 0009, 
"...estimated multi-fluid flow rates are provided for a plurality of selected well locations"), 
wherein the processor system is configured to allocate the flow rates bv inverting the 
sensed temperature data using a temperature forward model (paragraph 0039, lines 4- 
5; paragraph 0031 , lines 1-4, "The data path 31 supplies transient data to the model, 30, 
such as... temperature data 43 measured at multiple downhole locations"; paragraph 
0036, "Transient well data is measured... including pressure and temperature data in the 
wellbore 12 above each flow entry being produced. ..measurement above each flow 
entry is not required for the solution of the inverse problem... In step 106, the 
mathematical model of the wellbore 12 is run to calculate the expected pressure and 
temperature values at the downhole sensor locations, the expected volumetric phase 
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flow rates at the wellhead 22, and sensitivity coefficients of the model response to each 
phase flow rate at each fluid entry location"; paragraph 0037, lines 4-10, "...the model- 
calculated expected pressure and temperature values for various well locations of step 
106 are preferably compared with the measured temperature and pressure values 
obtained in step 104 with respect to those same well locations..."; paragraph 0038, 
"...the deviation between the calculated and measured quantities (of step 108) may be 
used with the sensitivity coefficients of the model (from step 106) to determine changes 
necessary. ..In this way, the modeling comparisons may be reiterated... until an 
approximate match (within acceptable tolerances) is obtained between the calculated 
well properties and related flow rates and the measured well properties and flow 
rates. ..if the measured volumetric phase rates and pressure and temperature readings 
are in tolerable agreement with the expected values predicted by the model 30, the final 
estimates of the multiphase flow rates are provided as shown in step 118"; claim 27, 
"...the multiphase fluid flow rate program further comprises a model inversion 
algorithm"; paragraph 0009, "...estimated multi-fluid flow rates are provided for a 
plurality of selected well locations"), wherein the inverting comprises using an 
optimization algorithm that solves an optimization problem for calculating the flow rates 
(Figure 3, elements 108, 110; paragraph 0006, lines 5-8, "The multi-phase fluid flow 
rates are estimated by iteratively comparing measured static and transient well 
conditions with the model for the well"; paragraphs 0037-0038, "...In step 108, the 
model-calculated expected pressure and temperature values for various well locations 
of step 106 are preferably compared with the measured temperature and pressure 
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values obtained in step 104 with respect to those same well locations. Thus at step 108, 
the actual transient data is compared with the calculated expectations of the model...", 
"...the deviation between the calculated and measured quantities (of step 108) may be 
used with the sensitivity coefficients of the model (from step 106) to determine changes 
necessary. ..In this way, the modeling comparisons may be reiterated... until an 
approximate match (within acceptable tolerances) is obtained between the calculated 
well properties and related flow rates and the measured well properties and flow 
rates. ..if the measured volumetric phase rates and pressure and temperature readings 
are in tolerable agreement with the expected values predicted by the model 30, the final 
estimates of the multiphase flow rates are provided as shown in step 118"), where the 
optimization problem minimizes an error between the sensed temperature data and 
corresponding calculated temperature data from the model (Figure 3, elements 108, 
112; paragraphs 0037 and 0038, "...the deviation between the calculated and measured 
quantities (of step 108) may be used with the sensitivity coefficients of the model (from 
step 106) to determine changes necessary for the estimate of phase flow rates at each 
well entry point. ..the modeling comparisons may be reiterated. ..until an approximate 
match (within acceptable tolerances) is obtained between the calculated well properties 
and related flow rates and the measured well properties and flow rates..."). 
18. As to Claim 20, Shah et al teaches: wherein the temperature forward model 
specifies temperature as a function of zonal flow rates {model in which temperature is a 
function of zonal flow rates: paragraph 0025, lines 4-6, "It is well known that each of the 
entries 24 has its own fluid phase (whether it be oil, water or gas), flow rate, 
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temperature and hydrocarbon mixture composition"; paragraph 0027, lines 1-14, "...the 
model 30 takes into account the conservation of energy and mass, and consequently 
simulates the evolution of the temperature of the flowing fluid... It is known that pressure 
and temperature of the flowing fluids change as they travel up or down a flow 
path... Therefore, model 30 should take such factors into consideration"; paragraph 
0029, lines 8-11, "The user of the invention may specify... the temperature at each entry 
point 24 within the well 10"; paragraph 0036, lines 9-12, "In step 106, the mathematical 
model for the wellbore 12 is run to calculate the expected pressure and temperature 
values at the downhole sensor locations..."). 

1 9. As to Claim 21 , Shah et al teaches: wherein the temperature sensor comprises a 
distributed temperature sensor (paragraph 0025, lines 16-24, "...a plurality of sensors 
27 are deployed to take measurements at the various production zones 26 or other 
points of interest inside the wellbore 12. The sensors 27 are preferably downhole 
temperature and pressure transducers coupled to computer 32. ..may include fiber optic 
distributed temperature sensing ("DTS") systems..."). 

20. As to Claim 22, Shah et al teaches: wherein the processor system is configured 
to allocate flow rates in a multi-layer, multi-phase liquid producing well (Figure 1; 
paragraph 0025, lines 10-13; paragraph 0026, lines 1-7). 

21 . As to Claim 25, Shah et al teaches: wherein the wellbore is oriented generally 
vertically (Figure 1; paragraph 0024, lines 11-15). 

22. As to Claim 38, Shah et al teaches: a sensor to measure a total flow rate of the 
wellbore at a wellhead (paragraph 0011, lines 7-9, "...measuring the transient 
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temperature, pressure and wellhead flow rate"; paragraph 0036, lines 8-9, "Volumetric 
flow rate measurements for each phase at the wellhead 22 are obtained..."; paragraph 
0002, lines 5-16, wherein it is known in the art to use these techniques to measure 
wellhead flow rate); wherein the processor system is configured to allocate, using the 
model, the total flow rate among the plurality of wellbore zones based on the sensed 
temperature data to allocate flow rates (paragraph 0031 , lines 1-4, "The data path 31 
supplies transient data to the model, 30, such as... temperature data 43 measured at 
multiple downhole locations"; paragraph 0036, "Transient well data is 
measured. ..including pressure and temperature data in the wellbore 12 above each flow 
entry being produced... measurement above each flow entry is not required for the 
solution of the inverse problem... In step 106, the mathematical model of the wellbore 12 
is run to calculate the expected pressure and temperature values at the downhole 
sensor locations, the expected volumetric phase flow rates at the wellhead 22, and 
sensitivity coefficients of the model response to each phase flow rate at each fluid entry 
location"; paragraph 0037, lines 4-10, "...the model-calculated expected pressure and 
temperature values for various well locations of step 106 are preferably compared with 
the measured temperature and pressure values obtained in step 104 with respect to 
those same well locations..."; paragraph 0038, "...the deviation between the calculated 
and measured quantities (of step 108) may be used with the sensitivity coefficients of 
the model (from step 106) to determine changes necessary... In this way, the modeling 
comparisons may be reiterated... until an approximate match (within acceptable 
tolerances) is obtained between the calculated well properties and related flow rates 
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and the measured well properties and flow rates. ..if the measured volumetric phase 
rates and pressure and temperature readings are in tolerable agreement with the 
expected values predicted by the model 30, the final estimates of the multiphase flow 
rates are provided as shown in step 118"; claim 27, "...the multiphase fluid flow rate 
program further comprises a model inversion algorithm"; paragraph 0009, "...estimated 
multi-fluid flow rates are provided for a plurality of selected well locations"). 

23. As to Claim 26, Shah et al teaches: a method, comprising: deploying a 
distributed temperature sensor along a wellbore (paragraph 0025, lines 16-24, "...a 
plurality of sensors 27 are deployed to take measurements at the various production 
zones 26 or other points of interest inside the wellbore 12. The sensors 27 are 
preferably downhole temperature and pressure transducers coupled to computer 
32. ..may include fiber optic distributed temperature sensing ("DTS") systems..."); 
utilizing a model of temperature as a function of fluid flow rates in the wellbore 
(paragraph 0025, lines 4-6, "It is well known that each of the entries 24 has its own fluid 
phase (whether it be oil, water or gas), flow rate, temperature and hydrocarbon mixture 
composition"; paragraph 0027, lines 1-14, "...the model 30 takes into account the 
conservation of energy and mass, and consequently simulates the evolution of the 
temperature of the flowing fluid... It is known that pressure and temperature of the 
flowing fluids change as they travel up or down a flow path... Therefore, model 30 
should take such factors into consideration"; paragraph 0029, lines 8-1 1 , "The user of 
the invention may specify... the temperature at each entry point 24 within the well 10"; 
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paragraph 0036, lines 9-12, "In step 106, the mathematical model for the wellbore 12 is 
run to calculate the expected pressure and temperature values at the downhole sensor 
locations..."); obtaining measured temperatures from the distributed temperature sensor 
(paragraph 0025, lines 16-24, "...a plurality of sensors 27 are deployed to take 
measurements at the various production zones 26 or other points of interest inside the 
wellbore 12. The sensors 27 are preferably downhole temperature and pressure 
transducers coupled to computer 32..."; paragraph 0031, lines 1-4, "The data path 31 
supplies transient data to the model 30, such as... temperature data 43... "Figure 3, 
element 104); determining fluid flow rates in corresponding wellbore zones using the 
measured temperatures in conjunction with the model (paragraphs 0036-0038; 
paragraph 0009, "...estimated multi-fluid flow rates are provided for a plurality of 
selected well locations"); wherein the determined fluid flow rates are calculated using an 
optimization algorithm that solves an optimization problem (Figure 3, elements 108, 
110; paragraph 0006, lines 5-8, "The multi-phase fluid flow rates are estimated by 
iteratively comparing measured static and transient well conditions with the model for 
the well"; paragraphs 0037-0038, "...In step 108, the model-calculated expected 
pressure and temperature values for various well locations of step 106 are preferably 
compared with the measured temperature and pressure values obtained in step 104 
with respect to those same well locations. Thus at step 108, the actual transient data is 
compared with the calculated expectations of the model...", "...the deviation between 
the calculated and measured quantities (of step 108) may be used with the sensitivity 
coefficients of the model (from step 106) to determine changes necessary... In this way, 
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the modeling comparisons may be reiterated... until an approximate match (within 
acceptable tolerances) is obtained between the calculated well properties and related 
flow rates and the measured well properties and flow rates. ..if the measured volumetric 
phase rates and pressure and temperature readings are in tolerable agreement with the 
expected values predicted by the model 30, the final estimates of the multiphase flow 
rates are provided as shown in step 118"), where the optimization problem minimizes an 
error between the measured temperatures and corresponding temperatures calculated 
by the model (Figure 3, elements 108, 112; paragraphs 0037 and 0038, "...the deviation 
between the calculated and measured quantities (of step 108) may be used with the 
sensitivity coefficients of the model (from step 106) to determine changes necessary for 
the estimate of phase flow rates at each well entry point.. .the modeling comparisons 
may be reiterated... until an approximate match (within acceptable tolerances) is 
obtained between the calculated well properties and related flow rates and the 
measured well properties and flow rates..."). 

24. As to Claim 27, Shah et al teaches: wherein determining the fluid flow rates 
comprises inverting the measured temperatures using the model to obtain the fluid flow 
rates (paragraph 0031 , lines 1-4, "The data path 31 supplies transient data to the 
model, 30, such as. ..temperature data 43 measured at multiple downhole locations"; 
paragraph 0036, "Transient well data is measured. ..including pressure and temperature 
data in the wellbore 12 above each flow entry being produced... measurement above 
each flow entry is not required for the solution of the inverse problem... In step 106, the 
mathematical model of the wellbore 12 is run to calculate the expected pressure and 
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temperature values at the downhole sensor locations..."; paragraph 0037, lines 4-10, 
"...the model-calculated expected pressure and temperature values for various well 
locations of step 106 are preferably compared with the measured temperature and 
pressure values obtained in step 104 with respect to those same well locations..."; 
paragraph 0038, "...the deviation between the calculated and measured quantities (of 
step 1 08) may be used with the sensitivity coefficients of the model (from step 1 06) to 
determine changes necessary... In this way, the modeling comparisons may be 
reiterated... until an approximate match (within acceptable tolerances) is obtained 
between the calculated well properties and related flow rates and the measured well 
properties and flow rates. ..if the measured volumetric phase rates and pressure and 
temperature readings are in tolerable agreement with the expected values predicted by 
the model 30, the final estimates of the multiphase flow rates are provided as shown in 
step 118"; claim 27, "...the multiphase fluid flow rate program further comprises a model 
inversion algorithm"). 

25. As to Claim 28, Shah et al teaches: wherein deploying the distributed 
temperature sensor comprises deploying the distributed temperature sensor in a 
generally vertical wellbore (Figure 1 , element 27; paragraph 0024, lines 1 1 -1 5). 

26. As to Claim 29, Shah et al teaches: wherein deploying the distributed 
temperature sensor comprises deploying the distributed temperature sensor in a 
deviated wellbore (Figure 1, element 27; paragraph 0024, lines 11-15). 
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27. As to Claim 31, Shah et al teaches: wherein determining the fluid flow rates 
comprises determining flow rates for a single-phase liquid producing well (paragraph 
0026, lines 1-7; paragraph 0025, line 5). 

28. As to Claim 32, Shah et al teaches: wherein determining the fluid flow rates 
comprises determining flow rates for a multi-phase liquid producing well (paragraph 
0026, lines 1-7; paragraph 0025, line 5). 

Claim Rejections - 35 USC § 103 

29. The following is a quotation of 35 U.S.C. 103(a) which forms the basis for all 
obviousness rejections set forth in this Office action: 

(a) A patent may not be obtained though the invention is not identically disclosed or described as set 
forth in section 102 of this title, if the differences between the subject matter sought to be patented and 
the prior art are such that the subject matter as a whole would have been obvious at the time the 
invention was made to a person having ordinary skill in the art to which said subject matter pertains. 
Patentability shall not be negatived by the manner in which the invention was made. 

This application currently names joint inventors. In considering patentability of 
the claims under 35 U.S.C. 103(a), the examiner presumes that the subject matter of 
the various claims was commonly owned at the time any inventions covered therein 
were made absent any evidence to the contrary. Applicant is advised of the obligation 
under 37 CFR 1 .56 to point out the inventor and invention dates of each claim that was 
not commonly owned at the time a later invention was made in order for the examiner to 
consider the applicability of 35 U.S.C. 103(c) and potential 35 U.S.C. 102(e), (f) or (g) 
prior art under 35 U.S.C. 103(a). 

30. Claims 8-10 and 33 are rejected under 35 U.S.C. 103(a) as being unpatentable 
over Shah et al as applied to claims 1 , 7 and 26 above, in view of Finsterle ("iTough2 
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User's Guide", Earth Sciences Division, Lawrence Berkley National Laboratory, 
University of California, May 2000). 

31 . Shah et al teaches a method of determining production rates in a well comprising 
inverting measured temperatures by applying a model of temperature as a function of 
zonal flow rates in a well to determine an allocation of production rates from different 
producing zones in the well, wherein inverting comprises determining a degree of 
certainty in the production rates allocated, allocating flow rate in at least one wellbore 
zone using temperature data in conjunction of the model, and determining error in the 
fluid flow rate and compensating for model error when inverting using the model to 
determine the fluid flow rates (Figure 3, elements 108, 112; paragraph 0038). 

32. Shah et al does not expressly teach: (claim 8) wherein determining the degree of 
certainty comprises determining a degree of error in the model, the method further 
comprising compensating for the determined degree of error in the model in performing 
the inverting; (claim 9) wherein determining a degree of certainty comprises determining 
a degree of error in the measured temperatures, the method further comprising 
compensating for the determined degree of error in the model in performing the 
inverting; (claim 10) wherein determining the degree of certainty comprises determining 
a degree of error in well parameter values, the method further comprising compensating 
for the determined degree of error in the model in performing the inverting; (claim 33) 
determining a model error, a measurement error, and a well parameter error; and 
compensating for model error, measurement error and well parameter error when 
inverting using the model to determine the fluid flow rates. 
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33. Finsterle teaches that in predicting multiphase fluid and heat flow in the 
subsurface by means of numerical simulation, errors in the conceptual model usually 
have the largest impact on model predictions, and assigning parameter values to the 
numerical model is likely to be tedious and challenging (page 2, paragraphs 1 and 2) 
and therefore teaches the iTOUGH2 computer program that provides inverse modeling 
capabilities for the TOUGH2 simulator (a numerical simulator for multidimensional, 
nonisothermal flows of multiphase, multicomponent fluids in porous and fractured 
media) that not only estimates model-related parameters by automatically calibrating 
TOUGH2 models to laboratory or field data, the information obtained by evaluating the 
sensitivity of the calculated system response with respect to certain input parameters 
can be used to study the appropriateness of a proposed experimental design and to 
analyze the uncertainty of model predictions (page 4, paragraph 1 ). The iTOUGH2 
program taught by Finsterle includes (claim 8, 33) determining a degree of error in the 
model (page 4, section (2), "error analysis is performed"; page 7, item (4), "...Model 
output and measured data are compared only at discrete points in space and time, the 
so-called calibration points", equation 1.5.3 and description); (claim 9, 33) determining 
a degree of error in the measurements (page 3, paragraph 2, lines 4-5; page 10, item 5; 
page 27, element 2.5.3.1); (claim 10, 33) determining a degree of error in model 
parameter values (page 4, section (2), "error analysis is performed"; page 8, items 7 
and 8, "objective function"); (claims 8, 9, 10, 33) and compensating for model error, 
measurement error and model parameter error when inverting (page 6, items 8-9; page 
8, items (8 and 9), "...find the minimum objecting function by iteratively updating the 



Application/Control Number: 10/575,030 Page 26 

Art Unit: 2123 

model parameters. Since the model output z(p) depends on the parameters to be 
estimated, the fit can be improved by changing the elements of parameter vector p..."; 
page 10, item 5 and page 1 1 , equation 1 .6.2 that shows the objective function to be 
minimized, S, includes measurement errors). 

34. Shah et al and Finsterle are analogous art since they are both directed to 
modeling multiphase flows in a subsurface by means of inverse modeling. 

35. It would have been obvious to one of ordinary skill in the art at the time the 
invention was made to modify the method of determining production rates in a well 
comprising inverting measured temperatures by applying a model of temperature as a 
function of zonal flow rates in a well to determine an allocation of production rates from 
different producing zones in the well, wherein inverting comprises determining a degree 
of certainty in the production rates allocated, allocating flow rate in at least one wellbore 
zone using temperature data in conjunction of the model, determining error in the fluid 
flow rate and compensating for model error when inverting using the model to determine 
the fluid flow rates as taught by Shah et al to further include (claims 8, 33) wherein 
determining the degree of certainty comprises determining a degree of error in the 
model; (claims 9, 33) wherein determining a degree of certainty comprises determining 

a degree of error in the measurements (measured temperatures); (claims 10, 33) 
wherein determining the degree of certainty comprises determining a degree of error in 
model (well) parameter values; (claims 8, 9, 10, 33) wherein determining error in the 
fluid flow rate comprises compensating for model error, measurement error and well 
parameter error as taught by Finsterle since Finsterle teaches that in predicting 
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multiphase fluid and heat flow in the subsurface by means of numerical simulation, 
errors in the conceptual model usually have the largest impact on model predictions, 
and assigning parameter values to the numerical model is likely to be tedious and 
challenging (page 2, paragraphs 1 and 2) and therefore teaches the iTOUGH2 
computer program that provides inverse modeling capabilities for the TOUGH2 
simulator (a numerical simulator for multidimensional, nonisothermal flows of 
multiphase, multicomponent fluids in porous and fractured media) that not only 
estimates model-related parameters by automatically calibrating TOUGH2 models to 
laboratory or field data, the information obtained by evaluating the sensitivity of the 
calculated system response with respect to certain input parameters can be used to 
study the appropriateness of a proposed experimental design and to analyze the 
uncertainty of model predictions (page 4, paragraph 1). 

36. Claim 1 1 is rejected under 35 U.S.C. 1 03(a) as being unpatentable over Shah et 
al as applied to claim 1 above, in view of Akin ("Analysis of Tracer Tests with Simple 
Spreadsheet Models", Computers and Geosciences, 27, pages 171-178, 2001). 

37. Shah et al teaches a method of determining production rates in a well comprising 
inverting measured temperatures by applying a model of temperature as a function of 
zonal flow rates in a well to determine an allocation of production rates from different 
producing zones in the well, wherein the inverting comprises an optimization algorithm 
(Figure 3 and description, paragraphs 0033, 0036-0038). 
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38. Shah does not expressly teach wherein using the optimization algorithm 
comprises utilizing a generalized reduced gradient optimization algorithm. 

39. Akin teaches a method of matching field data to predictions from computer 
simulation programs in tracer studies used for reservoir characterization that uses 
function evaluations rather than full simulator runs that results in a large reduction in 
computing time (Abstract, lines 4-7; section 1, paragraph 1, sentence 1; section 1, 
paragraph 2, lines 1-6 and lines 20-28), wherein after the reservoir models are 
implemented in a spreadsheet, the models were then matched to experimental data 
using the Generalized Reduced Gradient nonlinear optimization code to minimize the 
objective function (page 174, column 2, paragraph 4). 

40. Shah et al and Akin are analogous art since they are both directed to modeling a 
reservoir using inverse modeling techniques that match the models to experimental data 
using optimization techniques that minimize an objective function. 

41 . It would have been obvious to one of ordinary skill in the art at the time the 
invention was made to modify the optimization algorithm as taught by Shah to utilize a 
generalized reduced gradient optimization algorithm as taught by Akin since Akin 
teaches a method of matching field data to predictions from computer simulation 
programs in tracer studies used for reservoir characterization that uses function 
evaluations rather than full simulator runs that results in a large reduction in computing 
time (Abstract, lines 4-7; section 1, paragraph 1, sentence 1; section 1, paragraph 2, 
lines 1-6 and lines 20-28). 
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42. Claim 16 is rejected under 35 U.S.C. 103(a) as being unpatentable over Shah as 
applied to claim 12 above, and further in view of Curtis (US Patent 3,913,398). 

43. Shah et al teaches a method of determining flow rates in a well comprising 
determining flow rates of a plurality of liquid phases through each of a plurality of well 
zones via measured temperatures. 

44. Shah et al does not expressly teach wherein determining the flow rates 
comprises determining flow rates of fluid injected into each of the plurality of well zones. 

45. Curtis teaches that temperature data has been useful in the studies of secondary 
recovery of crude petroleum and that in the secondary recovery of crude petroleum, it is 
important not only to locate each permeable formation accepting fluid but also to 
determine the flow rates at which the fluid enters each formation (column 1 , lines 24-33) 
and teaches that determining flow rates of fluid injected into a plurality of well zones 
using temperature data is well known in the art (column 1, lines 34-50; column 2, lines 
27-41). 

46. Shah et al and Curtis are analogous art since they are both directed to 
determining flow rates of fluids (such as oil) produced in wells at each of a plurality of 
possible entry points through the use of temperature data. 

47. It would have been obvious to one of ordinary skill in the art at the time the 
invention was made to modify the method of determining flow rates in a well comprising 
determining flow rates of a plurality of liquid phases through each of a plurality of well 
zones via measured temperatures as taught by Shah to further include determining flow 
rates of fluid injected into each of the plurality of well zones as taught by Curtis since 
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Curtis teaches that temperature data has been useful in the studies of secondary 
recovery of crude petroleum and that in the secondary recovery of crude petroleum, it is 
important not only to locate each permeable formation accepting fluid but also to 
determine the flow rates at which the fluid enters each formation (column 1 , lines 26-33) 
and teaches that determining flow rates of fluid injected into a plurality of well zones 
using temperature data is well known in the art (column 1, lines 34-50; column 2, lines 
27-41). 

48. Claims 23 and 24 are rejected under 35 U.S.C. 103(a) as being unpatentable 
over Shah as applied to claim 1 9 above, in view of Tubel (US Patent 6,082,454). 

49. Shah et al teaches a system comprising a production completion in a wellbore 
(Figure 1, elements 14, 16, 18, 20, 22, 24, 28). 

50. Shah does not expressly teach that the production completion (claim 23) 
comprises an electric submersible pumping system or (claim 24) a gas lift system. 

51 . Tubel teaches spooled coiled tubing strings (production/completion strings) 
which include desired devices and sensors that are assembled and tested at the factory 
prior to the deployment of the string such that it substantially increases the quality and 
reliability of such strings and reduces the deployment or retrieval time (column 2, lines 
7-21), wherein a production completion according to the invention comprises devices 
typically used with a production completion such as an electrical submersible pump 
and/or a gas lift device according to the particular application (column 5, lines 29-36). 
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52. Shah and Tubel are analogous art since they are both directed to production 
completions deployed in a wellbore. 

53. It would have been obvious to one of ordinary skill in the art at the time the 
invention was made to modify the production completion as taught by Shah to further 
include an electric submersible pumping system and a gas lift system as taught by 
Tubel since Tubel teaches spooled coiled tubing strings (production/completion strings) 
which include desired devices and sensors that are assembled and tested at the factory 
prior to the deployment of the string such that it substantially increases the quality and 
reliability of such strings and reduces the deployment or retrieval time (column 2, lines 
7-21), wherein the devices include an electrical submersible pump and/or a gas lift 
device among others that are used with the production string according to the particular 
application (column 5, lines 29-36). 

Response to Arguments 

54. Applicant's arguments filed 6/10/10 have been fully considered but they are not 
persuasive. 

55. As to Claims 1,19 and 26, Applicant argues that Shah fails to disclose using an 
optimization algorithm that solves an optimization problem for calculating production 
rates from different producing zones in a well, where the optimization problem 
minimizes an error between the measured temperatures and corresponding 
temperatures calculated by the model. More specifically, Applicant argues that Shah's 
teaching of approximately matching calculated and measured quantities to within 
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acceptable tolerances does not provide any hint of an optimization algorithm that solves 
an optimization problem that minimizes an error between measured temperatures and 
corresponding calculated temperatures from a model (page 8). 

The Examiner respectfully disagrees. These arguments refer to limitations that 
have been amended into the claim. Accordingly, the Examiner has cited portions of the 
Shah reference for these limitations above. Shah's Figure 3, specifically, elements 108 
and 1 12, as well as paragraphs 0037 and 0038 clearly teach an optimization algorithm 
that minimizes the error between the measured and calculated temperatures. 
Specifically, element 108 compares the calculated and measured temperatures. This 
comparison produces a "deviation", that is, the "error" or "difference", between the 
calculated or measured temperatures. If the "deviation", or "error" is too large (does not 
fall within an acceptable tolerance level), changes are made to the model and the 
process is reiterated until the calculated and measured temperature deviation is within 
an acceptable tolerance level. These changes to the model and the reiteration of the 
process shown in Figure 3 and described in paragraphs 0037 and 0038 is clearly an 
optimization algorithm that is minimizing the error between the model calculations and 
actual measured data. 

56. As to Claim 1 2, Applicant argues that there is no teaching or hint in Shah of using 
a total flow rate from the well to allocate the total flow rate among the plurality of well 
zones in performing inverting to determine flow rates of a plurality of liquid phases 
through the plurality of well zones (page 9). 

The Examiner respectfully disagrees. These arguments refer to limitations that 
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have been amended into the claim. Accordingly, the Examiner has cited portions of the 
Shah reference for these limitations above. Shah teaches that the wellhead flow rate is 
measured and input to the model (paragraph 001 1 , lines 7-9; Figure 3, element 104; 
Figure 4, element 70). This measurement is input to the model and used to allocate the 
total flow rates of a plurality of liquid phases among a plurality of well zones (Figure 3; 
paragraphs 0036-0038). It is the Examiner's position that these cited portions of Shah 
clearly teach the limitations of the claim. 

Conclusion 

57. THIS ACTION IS MADE FINAL. Applicant is reminded of the extension of time 
policy as set forth in 37 CFR 1 .136(a). 

A shortened statutory period for reply to this final action is set to expire THREE 
MONTHS from the mailing date of this action. In the event a first reply is filed within 
TWO MONTHS of the mailing date of this final action and the advisory action is not 
mailed until after the end of the THREE-MONTH shortened statutory period, then the 
shortened statutory period will expire on the date the advisory action is mailed, and any 
extension fee pursuant to 37 CFR 1 .136(a) will be calculated from the mailing date of 
the advisory action. In no event, however, will the statutory period for reply expire later 
than SIX MONTHS from the mailing date of this final action. 

58. Any inquiry concerning this communication or earlier communications from the examiner 
should be directed to Mary C. Jacob whose telephone number is 571-272-6249. The examiner 
can normally be reached on Tuesday-Thursday, 7AM-4PM. 
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If attempts to reach the examiner by telephone are unsuccessful, the examiner's 
supervisor, Paul Rodriguez can be reached on 571-272-3753. The fax phone number for the 
organization where this application or proceeding is assigned is 571-273-8300. 

Information regarding the status of an application may be obtained from the Patent 
Application Information Retrieval (PAIR) system. Status information for published applications 
may be obtained from either Private PAIR or Public PAIR. Status information for unpublished 
applications is available through Private PAIR only. For more information about the PAIR 
system, see http://pair-direct.uspto.gov. Should you have questions on access to the Private 
PAIR system, contact the Electronic Business Center (EBC) at 866-217-9197 (toll-free). 



/Mary C Jacob/ 
Examiner, Art Unit 2123 
7/6/10 



